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Abstract. We develop a theory which describes hybrid structures consisting out 
of superconducting and ferromagnetic parts. We give two examples for applications. 
First, we consider a hybrid structure containing a strong ferromagnet in the ballistic 
limit. Second, we study for a weak ferromagnet the influence of a domain waU on 
the superconducting proximity effect. In both cases we account quantitatively for 
the mixing between singlet and triplet correlations. 

1 Introduction 

Hybrid structures containing ferromagnetic materials became recently a focus 
of nanoelectronic research because of their relevance for spintronics applica- 
tions. Consequently, it is desirable to understand how in the case of a super- 
conductor coupled to a ferromagnetic material superconducting correlations 
penetrate into the ferromagnet. A powerful method to treat such problems 
is the quasiclassical theory of superconductivity developed by Larkin and 
Ovchinnikov and by Eilenberger [1,2]. Within this theory the quasiparticle 
motion is treated on a classical level, whereas the particle-hole and the spin 
degree of freedoms are treated quantum mechanically. The paper consist of 
two parts. First we demonstrate our method for a ballistic heterostructure 
containing a strong ferromagnet, for which we chose for simplification a com- 
pletely polarized ferromagnet, a half metal. We show results for the modifi- 
cation of the quasiparticle density of states due to the proximity effect in the 
half- metallic material in a superconductor/half metal/superconductor device 
(see Fig. 1, left). Second, we study for a diffusive heterostructure containing 
a weak ferromagnet the influence of a domain wall of the Bloch type on the 
quasiparticle spectrum. The corresponding setup is shown on the right in Fig. 
1. We show results for different ratios between the domain wall width and 
the superconducting coherence length. 
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Fig. 1. Schematic picture of the studied devices. On the left, a ballistic heterostruc- 
ture where a half metal is sandwiched between two superconductors. Here, the con- 
ventional proximity effect is completely suppressed as a result of the complete spin 
polarization of the half metal. On the right, a diffusive weak ferromagnet between 
two superconductors. The ferromagnet contains a Bloch domain wall of size dw 



2 Basic Equations: Ballistic Case 



2.1 Transport Equations 

The central quantity in quasiclassical theory of superconductivity [1,2] is the 
quasiclassical Green's function .g(pj7, R, E, t) that depends on the spatial co- 
ordinate R and time t. It describes quasiparticles with energy E (measured 
from the chemical potential) and Fermi momentum pp moving along clas- 
sical trajectories with direction given by the Fermi velocity Vi?(pF).[3] The 
quasiclassical Green's function is a functional of self energies S{pf, R, E, t), 
which in general include molecular fields, the superconducting order param- 
eter Z\(pi7', R, t), impurity scattering, and external fields. The quantum me- 
chanical degrees of freedom of the quasiparticles show up in the matrix struc- 
ture of the quasiclassical propagator and the self energies. It is convenient 
to formulate the theory using 2x2 matrices in Keldysh space (denoted by 
a "check" accent), the elements of which in turn are 4x4 Nambu-Gor'kov 
matrices in combined particle-hole (denoted by a "hat" accent) and spin 
space. The structure of the propagators and self energies in Keldysh-space 
and particle-hole space is as follows. 



^ — { A V ^^'"^ - ( aR.a ^r,a ] > - ( aK ]■ (2) 
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The elements of the 2x2 Nambu-Gor'kov matrices are 2x2 matrices in 
spin space, e.g. = g^^^ with a,f3 = {T,i}, and similarly for others. In 
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writing Eqs. (1) and (2) we used general symmetries, which are accounted 
for by the "tilde" operation, 

X{pF,Ii,E,t)=X{-pF,Ii,-E,t)*. (3) 

The quasiclassical Green's functions satisfy the Eilenberger-Larkin-Ovchin- 
nikov transport equation and normalization condition 

[£f3-r,ff]^+jvF-Vff = 0, g^g = -TTn. (4) 

The noncommutativc product (x) combines matrix multipUcation with a con- 
volution over the internal variables, and fa = fal is a Pauli matrix in particle- 
hole space. 

The functional dependence of the quasiclassical propagator on the self 
energies is given in the form of self-consistency conditions. E.g for a weak- 
coupling, s-wave order parameter the condition reads 

A(R,t) = x ^(/^(pF,R,£;,t))p., (5) 

where A is the strength of the pairing interaction, and ( denotes averaging 
over the Fermi surface. The cut-off energy Ec is to be eliminated in favor of 
the transition temperature in the usual manner. 

When the quasiclassical Green's function has been determined, physical 
quantities of interest can be calculated. E.g. the equilibrium local density 
of states at position R for quasiparticles with momentum pp and energy E 
(measured from the Fermi level) reads 

iV(pf,R,£;) =iV^-i-Tr [f35^(p^,R,£;) -f35^(p^,R,ii;)] . (6) 

where Np is the density of states on the Fermi surface. 

For heterostructures, the above equations must still be supplemented with 
boundary conditions at the interfaces. As these conditions are non-trivial in 
quasiclassical theory, we present these conditions in the following chapter in 
detail. 

2.2 Boundary Conditions 

In order to formulate the boundary conditions at an interface between two 
materials, wc define first an auxiliary propagator on each side of the interface 
[4]. Then we relate this auxiliary propagator to the full propagator via a 
transfer matrix. This approach is completely equivalent to a full scattering 
matrix approach as we will show at the end of this paragraph. However, the 
repeated Andreev scattering processes across the interface complicate the full 
scattering matrix approach, and we found the current method much easier 
to implement and numerically more stable [5,6]. 
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We use for the auxiliary propagator the notation ^r"'" and g'"' , where the 
upper index a = {I, r} determines the side of the interface (left or right). The 
lower index denotes the direction of the Fermi velocity. Incoming momenta 
(index i) are those with a Fermi velocity pointing towards the interface, 
and outgoing momenta (index o) are those with a Fermi velocity pointing 
away from the interface. We formulate here the boundary conditions for clean 
surfaces, which conserve the parallel component of the Fermi momentum, 
P||. The auxiliary propagators as function of p||, energy and time t are 
solutions of the quasiclassical transport equation with the exact self energies, 
together with the normalization condition, however subject to the auxiliary 
boundary conditions, 

ff«'0(p||,i;,i) = 5S(p||) gr''{pii,E,t) 5«(p||), (7) 

where the following symmetries hold, 

^ro(Pii) = ^oi(Pii)^ = ^o"(Piir'- (8) 

These boundary conditions are formally equivalent to boundary conditions for 
an impenetrable interface, however with a surface scattering matrix S"- deter- 
mined from the reflection amplitudes of the full scattering matrix as explained 
below. Once these auxiliary propagators are obtained, the full propagators 
can be obtained directly, without further solving the transport equation, in 
the following way. We define hopping amplitudes f|o, T[g, f^l and T^j, for the 
four channels (left incoming to right outgoing, right incoming to left outgoing, 
left outgoing to right incoming, and right outgoing to left incoming) which 
conserve the momentum component parallel to the interface, py. With the 
help of these amplitudes we solve for the transfer matrices, t", for incoming 
trajectories from the following equations, 

if{p\\,E,t) = fr/(p||) gr{p\\,E,t) f„t(p||) 

® (l + gt'°{pll,E,t)^tf{pii,E,t)) , (9) 

where (a/3) = {{I, r), (r, /)}. The corresponding transfer matrices for outgoing 
trajectories are related to the ones for incoming trajectories through the 
relations 

i^{pil,E,t) = ^«(p||) it{pii,E,t) SZ{p\i). (10) 

Particle conservation requires certain symmetries between the hopping ele- 
ments, which are, 

^ot(Pll) = C^(Pll)^ = S^oM)rtiP\\)SoM)- (11) 

Consequently, only one amplitude, e.g. f|^, contains free material parameters, 
the other three amplitudes depend on it. The hopping amplitudes are defined 
via the transmission amplitudes of the full scattering matrix as is shown 
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below. Formally, they describe the modifications of the decoupled problem 
due to virtual hopping processes to the opposite side. 

The full propagators, fulfilling the desired boundary conditions at the 
interface, can now be easily calculated. For incoming trajectories they are 
obtained from 

+ (^gf'°{p\\,E,t)+iwl)(3tf{p\\,E,t)^(^gf'°{p\\,E,t)-inl), (12) 
and for outgoing trajectories from 

5„"(P||,i;,0=ffo"'°(P||,i^,t) 
+ {g^'\p\\,E,t) - iTvi) ^t'^ip\\,E,t) {g^'°{pi\,E,t) + inl) . (13) 

In this formulation, the boundary problem effectively reduces to calculating 
the auxiliary Green's functions for perfectly reflecting interfaces. Numerically 
this is an extremely simple task, e.g., employing the procedure of Riccati 
parameterization as explained below. Afterwards the boundary Green's func- 
tions for the partially transmitting interface can be obtained directly from 
Eqs. (12) and (13), since solving for the necessary transfer matrices (9) only 
involves a matrix inversion. 

In the transfer-matrix description, the phcnomcnological parameters con- 
taining the microscopic information of the interface are the two surface scat- 
tering matrices Sg[ and the hopping amplitude r-^. All three quantities are 
2x2 matrices in spin spac;e. As Sf^ is unitary, it depends apart from the direc- 
tion of the quantization axis on two parameters, a scalar scattering phase and 
a spin mixing angle (or, equivalently, a spin rotation angle). Similarly, as r-^ is 
hcrmitian, it also depends on two parameters, one describing spin conserving 
transmission and the other spin flip transmission. All remaining quantities 
defined above are related to these material parameters by symmetries. 

The particle-hole structures of the surface scattering matrix and the hop- 
ping amplitude are given by, 

SoM)~[ ) ' ^"^^^^ ~ \ fr/(P||)y'' ^''^ 

with the hole components, 

S:M) = 5f„(-p||)* = 5„",(-p||)*'^ (15) 

^ro''(Pii) = r:f{-p\0* = s^A-pwT rr/(-Pii)* ^f,(-Pii)*. (i6) 

Finally, we relate the parameters and T^f to the full normal state 
scattering matrix S, 
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The scattering matrix is diagonal in particle-hole space, with diagonal com- 
ponents 

S„„ = (1 + ^Xtri")-' (1 - nXtrf:) 5", (18) 

and off-diagonal components 

S^p = {l + ^Xtrt-)-'2nr:t. (19) 

These identities serve as a precise definition of the auxilary parameters of the 
theory, 5oj(p||), S^^{^p\\)^ and t]^(P||) in terms of the physical parameters of 
the full scattering matrix. 



2.3 Riccati parameterization 

The method of the Riccati parameterization [7-9] of the quasiclassical Green's 
functions has proved powerful during recent years. It accounts automatically 
for the normalization condition. A corresponding parameterization for distri- 
bution functions was applied to the Keldysh part of the Green's functions. 
The combined equations in their general form are [9] 

Y-(7-" (g) a; - a; 7^) (x - 7-" x (g) 7^) I ' 



g -±nrN _(i ^ ^ , (20) 

with 

-(^ (l-7«.^®7«.^)-iJ • ^^^^ 

Thus, the problem is reduced to the solution for two 2x2 matrices in spin 
space, 7^ and x. Using fundamental symmetries (particle-hole, retarded- 
advanced) like Eq. (3) and 

7^(PF, R, E, t) = 7''(-PF, R, -E, ^)*^ (22) 

where tr denotes the spin-matrix transpose operation, the full retarded, ad- 
vanced and Keldysh Green's functions are obtained from 7^ and x. The 
transport equations for the functions j^' {p p, II, E,t) and x{pp,Il, E,t) are 

2£;7^ + jvirV7^ = 7^ (g) (g) 7^ -h (g) 7^ - 7^ (g) - A^, (23) 

and 

idtx + iwp^x + (-7^ (g - r^) (g X + X (g {-A^ (g 7-^ + T^) 

= -7^(gr^(g7^-hZ\^(g7^+7^(g^^-i:^. (24) 
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3 Results: Ballistic Case 

As an example we consider here the Andreev quasiparticle spectrum in a 
half-metaUic ferromagnet between two singlet superconductors. This corre- 
sponds to the device shown on the left side in Fig. 1. In this case, on the 
half-metallic side only quasiparticles with one spin direction with respect to 
the quantization axis are itinerant. The stable equilibrium configuration of 
such a system is a 7r-junction, in which the phase of the singlet order pa- 
rameter on both sides of the interface differs by tt [5] . Characteristic of such 
a heterostructure is the presence of triplet correlations near the interface. 
Shown in Fig. 2 are the singlet order parameter and all three triplet compo- 
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Fig. 2. Self consistent order 
parameter and triplet cor- 
relations in an S/HM/S tt- 
junction. The calculations are 
for temperature T — 0.05Tc, 
and for 'rj^,t/''"T,T ~ ^■'^ 



nents which are induced near the interface by the spin rotation effect [10]. 
The spin rotation effect is taken into account via a surface scattering ma- 
trix Soi = exp{i6az/2)l at the superconducting side of the interface, where 
6 defines a spin-rotation angle and CTz denotes the Pauli spin matrix [10,11]. 
Generally, the value of depends on the angle of impact, ip [10] and can 
approach values of the order of tt for strong band splitting [12]. For definite- 
ness, we present results for 9 — 0.757r cos^". On the half-metallic side the 
scattering matrix has no spin structure. Further, we use the hopping ampli- 
tude Tio = {1 + Sio)To cos ip, where tq = ('''T.T' ''"ij)'^ determined by the two 
spin scattering channels from the superconductor to the half-metallic spin- 
up band. This reflects a spin rotation during transmission which is half of 
the spin rotation during reflection. The cos^ factor accounts for the reduced 
transmission at large impact angles. We present results for T|,t/'''T,T ~ 0.7 
and 0.1, 27rT|.| = 1.0. The length of the half metallic region is 3^o with the 
coherence length — w//27r Tc. We assume cylindrical Fermi surfaces. 
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We have studied the behavior of such a heterostructure as a Josephson 
device. If a phase difference cj) in the order parameter is present between the 
left and right end of the heterostructure, the quasiparticle spectrum in the 
half metal is modified. In Figs. 3 and 4 we show results of our calculations 
for two sets of parameters. 




-1 -0.5 0.5 1 -1 1 



Fig. 3. Density of states at T = 0.05Tc for quasiparticles with normal impact at 
the half-metallic side of the left interface, for r|,i-/'''T,T = O-^- a) Dispersion of the 
maxima of the density of states as function of phase difference. Gapped regions 
are white, ungapped regions gray. The signs indicate the direction of the current 
carried by the Andreev states, b) As an example we show spectra for a selected 
phase difference 

The spectra consist out of Andreev quasiparticle bands, separated by 
gaps. In Figs. 3 a) and 4 a) we show these gapped-regions (white) and the non- 
gapped regions (grey) of the density of states at the half-metallic side of the 
left interface as a function of the phase difference. One characteristic spectrum 
for fixed phase difference for both positive and negative momentum direction 
is shown in Figs. 3 b) and 4 b). The dispersion of the maxima in the density 
of states with phase difference determines the sign of the Josephson current 
trough the device. We indicate in Figs. 3a) and 4a) by -|- and — Andreev 
states, which carry current in positive and negative direction respectively. 
It can be seen that the sign of the Josephson current corresponds to the 
sign of dE{(p)/d(p, where E((p) denotes the dispersion of the maxima in the 
density of states. This is analogous to the current carried by a bound state 
Ef, dispersing with phase (j), which is given by Jb ~ {dEb/d(j))n(Eb), where 
n{Eb) is the equilibrium fermion distribution function. 



Singlet-Triplet Mixing in SF hybrid devices 



9 



0.5 



-0.5 



-1 



a) 



-1 -0.5 0.5 1 -1 

E/A„ E/A„ 




Fig. 4. The same as Fig. 3 for ti^t/t^,]- = 0.1 



4 Basic Equations: Diffusive Case 



4.1 Transport Equations 

The fundamental quantity for diffusive transport is the Usadel Green's fune- 
tion [13], which is the momentum average of the quasiclassical Green's func- 
tion g(R, E, t) = {g{pF, R- E, t))pp . It is a fimctional of momentum averaged 
self energies S(R., E, t) — (^^(pi?, R, E, t))pp. The structures of g and S are 
the same as in Eqs. (1) and (2). Eq. (3) is replaced by 



X{Ii, E,t) =X(R,-E,ty. 



(25) 



The Usadel Green's function obeys the following transport equation and nor- 
malization condition [13], 



[ET3i-E,g]^ 



D. 



-Vi {g ( 



i Vjg) = 0, 



(26) 



Here, summation over the repeated index j is implied. S does not contain 
the non-magnetic impurity scattering self energy anymore. We will solve these 
equations by using the Riccati parameterization technique and present below 
the resulting equations for the diffusive case. The corresponding equations in 
a finite vector potential are obtained by applying a gauge transformation to 
Eqs. (26) (and to Eqs. (28), (29) below). 



4.2 Boundary Conditions 

In the diffusive limit the boundary conditions are formulated in terms of the 
momentum averaged transfer matrices and Green's functions. The resulting 
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boundary conditions for spin- conserving interfaces were found by Nazarov 
[14] and follow from the transfer-matrix approach outlined above [6]. They 
are 



SRb 47r2-T({5/3,5«} + 27r2) 



az 



(27) 



where z is the coordinate along the interface normal, (7a(f3) and g"^^^ refers 
to the conductivity and the Keldysh Green's function on side a{(3) of the 
interface, S is the surface area of the contact, Rb the boundary resistance 
and T the transmission probability, T = p^^Jppp-- 



4.3 Riccati parameterization 

We use the same parameterization, Eqs. (20) (21), for the momentum av- 
eraged Green's functions. Note that the such defined diffusive quantities 
•y^CRjEjt) and a;(R, i^, t) are by no means in a simple way related to the 
ballistic quantities j^'{pf:'R., E, t) and x{pp,Ii, E, t). Consequently, also the 
structure of the transport equations is very different. In particular the gra- 
dient terms in the transport equations for these quantities are modified with 
respect to the gradient terms in the transport equations for the ballistic case. 
Instead of Eqs. (23) and (24) we obtain for the diffusive case 

fR 1 



2£;7" - iD 



+ (V,-7'') 



7 



^7^ 

)7«. 



■7 



(28) 



and 



idtx — iD 



f ®(V,-7^) 



+ (V,-7^) 
+ ( 



R 



-7 



7 



in 



{Vjx) + (Vjx) 
E^) X + X IS) {■ 



— 

ITT 



7^ + S^) 



7^ + Zi^ 



1 7^ + 7^ ' 



(29) 

Summation over the repeated index j is implied. The expressions for f^, 
and g^ are obtained by comparing Eq. (1) with Eqs. (20)~(21). 



5 Results: Diffusive Case 

As an example we study for this case the heterostructure shown on the right 
side of Fig. 1, which contains a weak diffusive ferromagnet between two su- 
perconductors. A domain wall of Bloch type is centered in the ferromagnetic 
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part. Problems involving domain walls have previously been treated only 
with techniques that do not allow access to information on spatial variations, 
[15] or elaborate methods involving rotating coordinate systems. [16] However, 
since the Riccati description contains the full 2x2 spin structure, it extends 
easily to account for any such situation. 

We adopt for the description of a weak ferromagnetic materials {J much 
smaller than the Fermi energy) a quasiparticle dispersion given by a spin- 
dependent energy shift Ef^ Ef^ — J(R) ■ a, in Eq. (26). Here J(R) 
denotes the effective exchange field of the ferromagnet, and cr denotes the 
vector of Pauli spin matrices. We solve the Usadel equations using the Riccati 
parameterization for a ferromagnetic material, which contains a domain wall 
of the Bloch type. We model the domain wall by a J- vector rotating across the 
domain wall with the polar angle 9j = a,ta,n{z / dw) , where dw parameterizes 
the width of the domain wall. 




|J|=0.0 

0.5 A, cl_W=0.2 5 
d_W=0.4 t, 
d_W=0.8 i 
d_W=1.0 5 
d_W=2.0 t. 



Fig. 5. Local density of states 
as a function of energy in 
the center of a heterostruc- 
ture containing a weak ferro- 
magnet (J = 0.5A) of length 
dp ~ 2^0 between two super- 
conductors. A Bloch domain 
wall of length dw is centered 
in the ferromagnet (see Fig. 
1). Results for various d^y are 
shown 



In Fig. 5 we present results for the local density of states at the center 
position of the heterostructure shown in Fig 1. The length of the ferromagnet 
is dp = 2^0: where = \/D/^ is the coherence length in the ferromagnet, 
and A is the gap in the superconductors. The dotted curve shows the be- 
havior for a normal metal between two superconductors, showing the well 
known minigap. The effect of a domain wall in the center of the ferromagnet 
is shown for several wall thicknesses dw in Fig. 5. The new structures in the 
presence of the domain wall in Fig. 5 are due to spin-triplet correlations which 
mix between all three spin-triplet channels due to the continuous variation 
of the spin quantization axis across the domain wall [16]. This leads to the 
appearence of additional Andreev bound states inside the minigap, modify- 
ing the total density of states. In general, the presence of the domain wall 
reduces the minigap, and for thick enough domain walls the minigap closes 
completely. This suggests that the width of the domain walls can influence 
the transport properties of such devices considerably. 
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6 Conclusions 

We have developed a framework for studying heterostructures with super- 
conducting and ferromagnetic parts, and have appUed it to two cases: a fer- 
romagnet in the baUistic regime with strong spin polarization (J ^ Ep): and 
a weak ferromagnet (J <ti Ep) in the diffusive regime, having an inhomoge- 
neous, non-collinear spin magnetization (Bloch domain wall). In both cases 
we have studied the superconducting proximity effect with the ferromagnet 
sandwiched between two superconductors. We have shown, that singlet-triplet 
mixing occurs near interfaces between superconductors and strong fcrromag- 
nets, or near domain-wall structures in the case of weak ferromagnets in 
proximity to a superconductor. The origin of equal-spin triplet correlations 
in the two studied cases is very different. In a strong ferromagnet between 
two superconductors, the singlet triplet mixing takes place within the super- 
conductor, in a layer of the order of the coherence length near the interface, 
enabling triplet correlations to penetrate into the ferromagnet. In contrast, 
in weak ferromagnets with domain walls of Bloch type, equal-spin triplet 
correlations are produced within the ferromagnet, near the domain walls. 
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